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ABSTRACT 

We report the detection of a 3.6a 350 /im-selected source in the Bootes Deep Field. The source, the first Short- 
wavelength Submillimeter-selected Galaxy (SSG 1), was discovered as part of a blank field extragalactic survey 
using the 350 /im-optimised Submillimeter High Angular Resolution Camera (SHARC II) at the Caltech Sub- 
milUmeter Observatory. With multiwavelength photometry from NOAO-NDWFS (R and I band), FLAMEX (J 
and Kj), Spitzer (IRAC and MIPS) and the Westerbork 1 .4 GHz Deep Survey (radio upper limit), we are able 
to constrain the photometric redshift using different methods, all of which suggest a redshift ^ L In the ab- 
sence of long-wavelength submillimeter data we use SED templates to infer that this source is an ultraluminous 
infrared galaxy (ULIRG) with a dust temperature of 30±5 K, occupying a region of luminosity-temperature 
space shared by modarate redshift /S'O-selected ULIRGs (rather than high redshift SCUBA-selected SMGs). 
SHARC II can thus select SMGs with moderately "warm" dust that might be missed in submillimeter surveys 
at longer wavelengths. 

Subject headings: infrared: galaxies - submillimeter: galaxies - galaxies: starburst - galaxies: high-redshift 



L introduction 

Submillimeter-selected galaxies (SMGs) generally refer to 
the population detected in the pioneering lensed and blank 
deep field surveys using the Submillimeter Common User 
Bolometer AiTay (SCUBA; Holland et al. 1999) instrument 
on the James Clerk Maxwell Telescope (JCMT) (e.g., Smail, 
Ivison & Blain 1997; Hughes et al. 1998; Barger et al. 
1998; Eales et al. 1999). They are regarded as the high red- 
shift (z-2-3; Chapman et al. 2005; Chapman et al. 2003; 
see also Simpson et al. 2004) counterparts to the popula- 
tion of luminous and ultraluminous infrared galaxies (LIRGs 
and ULIRGs) detected by the InfraRed Astronomical Satel- 
lite, IRAS, typically out to z<0. 1 (Soifer et al. 1 984; Joseph & 
Wright 1985; and Soifer, Neugebauer & Houck 1987), since 
they share similar properties (fundamentally, that the bulk of 
the bolometric luminosity is emitted in the restframe far-IR, 
powered by a combination of star formation and AGN activ- 
ity.). 

As a necessary consequence of the selection at longer sub- 
millimeter wavelengths, SMGs tend to have higher bolomet- 
ric luminosities, and those with comparable luminosities to 

' Imperial College London, Blackett Laboratoiy, Prince Consort Road, 
London SW7 2AZ, UK 

^ Observational Cosmology Laboratory (Code 665), NASA / Goddard 
Space Flight Center, Greenbelt, MD 20771 

3 SSAI, Lanham, MD 20706 
NRC 

' Steward Observatory, University of Arizona, 933 N. Cherry Avenue, Tuc- 
son, AZ 85721 

Spitzer Science Center, MC 220-6, 1200 East California Boulevard, 
Pasadena, CA91125 

' Department of Astronomy, Cornell University, 610 Space Sciences 
Building, Ithaca, NY 14853-6801 

* Harvard Smithsonian Center for Astrophysics, 60 Garden Street MS-66, 
Cambridge, MA 02138 

' National Optical Astronomy Observatory, Tucson, AZ 85719-6732 

'" Jet Propulsion Laboratory, California Institute of Technology. MC 169- 
327, 4800 Oak Grove Drive, Pasadena, CA 91 109 

' ' Department of Astronomy, University of Florida, Gainesville, FL 326 1 1 

'^ Deceased 



/RA5-selected ULIRGs tend to have cooler dust temperatures 
(Blain et al. 2004). Between the two are the Infrared Space 
Observatory (/50)-selected ULIRGs (Aussel et al. 1999; 
Rowan-Robinson et al. 1999; Elbaz et al. 2002a), which 
occupy an intermediate redshift space (z'^l), and have dust 
masses in-between the IRAS and SCUBA-selected ULIRGs 
(Chapman et al. 2002; Blain et al. 2004; this is also a re- 
gion of parameter space that will be shared by S/^/fzer-selected 
galaxies, see, e.g., Yan et al. 2004). 

Sources selected at 350 /im are expected to be predom- 
inandy LIRGs and ULIRGs at 1< z <3 (the K-correction 
changes from negative to positive with increasing redshift as 
the observing window moves towards and past the Wien peak; 
see, e.g., Guiderdoni et al. 1998; Khan, 2005, in preparation). 
Galaxies at z~l are likely to be the dominant source of the 
integrated far-IR background (Puget et al. 1996; Fixsen et al. 
1998; Elbaz et al. 2002b); also, observations suggest that after 
a relatively steep rise from z=0 to ~1, the cosmic star forma- 
tion rate (CSFR) appears to flatten off between z=l and ^A, 
with evidence for a slow decline at higher redshift (Lilly et al. 
1996; Madau et al. 1996; Connolly et al. 1997; Steidel et al. 
1999; Gabasch et al. 2004). Hence 350 /im-selected sources 
are an important probe of the epoch of peak star formation 
rate in the universe. 

In this Letter we report on the detection of a single source 
above 3ct at 350 /xm. It is, to the best of our knowledge, the 
first short-wavelength submillimeter-selected galaxy (200- 
500 /xm) and is denoted as SSG 1 in this Letter (but by stan- 
dai-d convention, SMM J143206.65H-341613.4). SSG 1 is the 
first object discovered purely by its 350 pm emission, through 
blank deep observations of the Bootes field in a program de- 
signed to search efficiently for new sources at 350 pm (the 
SHARC II Unbiased Deep Survey, SUDS; Khan et al., 2005, 
in preparation). Verification of the detection is confirmed 
by coincident sources subsequently identified in the Spitzer 
MIPS 24 /im filter, all four Spitzer IRAC filters, and J, K,, I 
and R band, enabling us to place constraints on the photomet- 
ric redshift, luminosity and nature of this source. This Letter 
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Table 1 

multiwavelength photometry for ssg 1 



Observed wavelength 


Flux density 


Instrument-Survey 


X-ray (0.5-7 keV) 


< 4 X 10^'^c.g.s 


C/ianrfra-XBootes 


4220A 


< 0.1 fiJy 


KPNO-NDWFS 


6590A 


1.0±0.1/xJy 


KPNO-NDWFS 


8081A 


3.1 ±0.2/xJy 


KPNO-NDWFS 


1.24 fim 


10.5 ± 1.1 /tJy 


KPNO-FLAMEX 


2.16 fim 


42.0 ± 1.4 ^ijy 


KPNO-FLAMEX 


3.6 


80.0 ± 2.9 AiJy 


Spitzer IRAC-Shallow 


4.5 fim 


65.8 ± 3.6 ^tJy 


Spitzer IRAC-Shallow 


5.8 


55.3 ± 17.3 AtJy 


Spitzer IRAC-Shallow 


8.0 Atm 


81.7 ± 15.1 /xJy 


Spitzer IRAC-Shallow 


24 /im 


0.61 ± 0.04 mJy 


Spitzer MIPS-IRS/MIPS GTO 


70 fj,m 


< 40 mJy 


Spitz.er MIPS-IRS/MIPS GTO 


160 /im 


< 100 mJy 


Spitzer MIPS-IRS/MIPS GTO 


350 /im 


23.2 ± 7.9 mJy 


SHARC II-SUDS 


21 cm 


< 63 /iJy 


Westerbork-Deep 1.4 GHz 



presents these fluxes, plus important flux upper limits in other 
bands, to characterize this source, and with this understand- 
ing we consider the implications for other sources that may 
be discovered at 350 /xm. 

2. OBSERVATION 

Eight hours of data were obtained using the Second 
Generation Submillimeter High Angular Resolution Cam- 
era (SHARC II) at the Caltech Submillimeter Observatory 
(CSO) on Mauna Kea, Hawai'i, in January and March 2004. 
SHARC II is a 350 /im-optimized camera (Dowell et al. 2003) 
built around a 12 x 32 element close-packed bolometer ar- 
ray (Moseley et al. 2004). It achieves a point-source sensi- 
tivity of ~ 1 Jy Hz"'/^ in good weather The 384 pixels of 
the SHARC II array image a region of around 1 .0' x 2.5' on 
the sky. Its filled absorber array provides instantaneous imag- 
ing of the entire field of view, sampled at roughly 2.5 pixels 
per nominal beam area. The beam profile was measured on 
known compact sources, and was verified to be within three 
per cent of the diffraction-limited beamwidth of 8.5". All ob- 
servations were taken using the Dish Surface Optimisation 
System (Leong et al. 2003), which corrects for the primary 
mirror deformation as a function of zenith angle, to improve 
the telescope efficiency and the pointing. 

For these data the in-band zenith atmospheric opacity 
(7350 /jm) ranged from 1.0 to 1.4, corresponding to a zenith 
transmission of around 30 per cent. Our observations were 
centred on the Bootes Deep Field (de Vries et al. 2002), at 
position RA= 14''32'"5.75^ Dec= 34°16'47.5" (J2000). 

The data were reduced using the standard CSO reduction 
software, CRUSH (Kovacs, 2005, in preparation). This soft- 
ware implements a self-consistent least-squares algorithm to 
solve for the celestial emission, taking into account instru- 
mental and atmospheric contributions to the signal. 

The skymap is calibrated with the flux and point spread 
function based on observations of Callisto taken throughout 
the observing period at similar elevations. 

3. RESULTS 

An oversampled At is used to determine the position of 
the source and the flux per beam. The best-fitting position is at 
14:32:06.65±0.24h-34:16:13.4±3.4 (J2000) (quoting 3 cr po- 



table 2 

Coincident detection positions and probability of 
accidental overlap within the 99% confidence region 

OF SSG 1 



Position (J2000) 


Band 


P(Accidental Overlap) 


RA Dec 






14:32:06.61 +34:16:12.0 


R 


9.5% 


14:32:06.55 +34:16:11.6 


K, 


2..9% 


14:32:06.58 +34:16:12.0 


3.6 iim 


9.4% 


14:32:06.60 +34:16:11.5 


24 pm 


0.2% 


14:32:06.65 +34:16:13.4 


350 fim 


n/a 



sition uncertainties), with a flux per beam of 23.2 ± 7.9 mJy 
(the 1(7 uncertainty comprises the background noise, primary 
calibrator and gain calibration errors). The signal-to-noise in 
the detection is 3.6. The additional photometry for SSG 1 
(Table 1) is from the Chandra XBootes survey (Murray et 
al. 2005), NOAO-NDWFS (Jannuzi & Dey 1999), FLAMEX 
(Elston et al., 2005, in preparation), Spitzer-IRAC Shallow 
survey (Eisenhardt et al. 2004), Spitzer-MlPS (IRS and MIPS 
GTO teams) and the Westerbork deep 1.4 GHz survey (de 
Vries et al. 2002). The upper limits are quoted to 3a. The 
corresponding multiwavelength positions and probability of 
coincident detection are given in Table 2. 




14''32"08'0 07"0 06'0 

RA (J3000) 



Fig. 1. — 350/.tm continuum emission map of SSG 1, as imaged by 
SHARC II. The SHARC II (plus) and MIPS position (cross) are overplot- 
ted, along with the resolved IRAC 3.6 pm positions of SSG 1 (open square) 
and SSG IE (filled square). The contours are in levels of 5mjy/beam. 

Within the 99% 350 /xm confidence region there are two 
candidate sources resolved in the optical and NIR bands (Fig- 
ure 1). We assume that the optical-NIR counterpart of SSG 1 
is the source coincident with the MIPS 24 /im detection. The 
source to the east of this (at IRAC position RA= 14''32"'6.94', 
Dec= +34°16'11.6" /J2000) is not detected by MIPS, and is 
designated as SSG IE for the purposes of this Letter. 

4. CONSTRAINING THE NATURE OF SSG 1 

Using the fluxes presented in Table 1, we can use several 
methods to constrain the redshift, reddening, infrared lumi- 
nosity and dust temperature of SSG 1. It can be seen in Fig- 
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Fig. 2.— Best-fitting observing frame SED to SSG 1 from STARDUST2 



Table 3 

Best-fitting PHOTOMETRIC redshifts and extinction for 
SSG 1 AND SSG IE WITH ImpZ and STARDUST2 



Source 


Photometric Redshift 


Model 


Av 


SSG 1 
SSG IE 


n QQ+0.04 

"•^^-0.03 

1.0±0.05 


IMPZ 
STARDUST2 
IMPZ 


2.6±0.2 
2.2±0.09 
0.5±0.1 



ure 2 that the optical-NIR SED of SSG 1 is characterized by 
a prominent bump associated with the continuum emission of 
the stellar populations. This bump peaks at 1 .6 /zm in the rest- 
frame, providing a very good constraint on the redshift, and 
also indicates SSG 1 is more likely to be starburst rather than 
AGN-dominated (which are usually associated with a feature- 
less power-law spectrum; see, e.g., Egami et al. 2004). 

ImpZ (Babbedge et al. 2004) uses only the optical-NIR pho- 
tometry when fitting templates with Bayesian statistics. The 
best-fitting redshifts are z=1.0![[;^5 SSG 1 and z=1.0±0.05 
for SSG IE respectively (Icr), and are listed in Table 3. It is 
worth noting that if, as ImpZ suggests, both galaxies are at the 
same distance, and interacting, it would make the SSG pair a 
widely separated interacting ULIRG system (a local example 
of which is IRAS 09111-1007; Khan et al. 2005). 

STARDUST2 (Chanial et al., 2005, in preparation) uses a 
mid-IR to radio spectral library constrained by the IR-radio 
correlation and a variety of local, IRAS, ISO and SCUBA 
color-color correlations in addition to a stellar synthesis 
model for the FUV to NIR window (Devriendt et al. 1999). 
This provides a simultaneous constraint on the thermal dust 
emission as well as the redshift. The effective dust tempera- 
ture is found following the methodology given in Chapman et 
al. (2002; 2005). The model assumes Hq = TOkms"' Mpc"', 
rim = 0.3, and ^Ia = 0.7, and that a single component is re- 
sponsible for both the optical and IR emission. Because of 



the degeneracy in fitting Tdust and the Total-IR Luminosity, 
Ltir (8-1 100 /im), it is the 350 /im flux in conjunction with 
the radio upper limit that constrains the infrared luminosity. 
The radio upper Umit (Table 1) is used to find as follows: 
X^ = J2xj with the radio xf term being 

xLo = [^(l•4GHz)-3^]VfT^/.(1.4GHz) > 3a 
= 0, /„(1.4GHz)<3<j 

where fi, is the modeled spectral energy distribution and 3a 
is the radio upper limit listed in Table 1 . 

The best-fitting model for SSG 1 from STARDUST2 re- 
turns a redshift of 0.99!jj;j]^ (all errors are quoted to 1 a; note: 
STARDUST2 uses the longer wavelengths, in addition to the 
optical-NIR photometry, to fit the extinction given in Table 3). 
This model also gives a dust temperature of 30.3±4.5K, and 
a log(LTiR) of 12.02![;;22Lq, with 88% of the total bolomet- 
ric luminosity radiated at wavelengths longward of 5 fim. The 
predicted 850 /im flux is 2.7^g ^ mJy. 

Other redshift estimates were independently obtained fol- 
lowing an approach that mostly reUes on the fit of the 1 .6 /xm 
feature (see, e.g., Le Floc'h et al. 2004). In agreement with 
the results derived from ImpZ and STARDUST2, these fits of 
the stellar bump led to a redshift z^l using the photometric 
Arp220 SED, and z= 1.25 ±0.25 using various templates from 
the library of Devriendt et al. (1999). 

We also use the radio- submiUimeter correlation to estimate 
a minimum redshift for the source, as the radio flux is un- 
likely to be significantly AGN-enhanced. Formally, these cor- 
relations should be used as a statistical redshift indicator, so 
with that caveat we present the minimum redshift of SSG 1 
based on the 1.4 GHz upper limit (Table 1) and an 850 /xm 
flux derived from flie best-fitting STARDUST2 SED (Figure 
2). Using the relations of Dunne, Clements & Eales (2000) 
and CariUi & Yun (2000a; 2000b) we get z„^n of 1.2 and 1.5 
respectively. 

The photometric redshift is more secure than the radio- 
submilhmeter correlation since the latter is prone to system- 
atics caused by uncertain dust temperatures (see Clements et 
al. 2004). 

5. discussion 

The (R-Ks)vega color of 5.7 for SSG 1 classifies it as an Ex- 
tremely Red Object (ERO), and bofli ImpZ and STARDUST2 
find this object to be highly reddened (Table 3). At least a 
third of the SMG population can be classified as EROs (Small 
et al. 2002; Webb et al. 2004; Frayer et al. 2004), which 
comprise two classes of galaxies: elliptical and dusty star- 
forming luminous infrared galaxies (but a significant submil- 
limeter flux is usually indicative of the latter class). EROs are 
thought to comprise a significant fraction of the cosmic star 
formation density at redshifts of one and higher (Cimatti et al. 
2002), the epoch by which the majority of the universe's star 
formation has taken place (Dickinson et al. 2003; Rudnick et 
al. 2003). 

Determining the effective dust temperature of the source al- 
lows us to estimate the total fraction of the bolometric lumi- 
nosity that is reprocessed by the dust. Our effective dust tem- 
perature from STARDUST2 is 30.3±4.5K, which, with our 
total infrared luminosity of 1 .0 x lO^^L©, allows us to directly 
compare SSG 1 with IRAS and SCUBA-selected ULIRGs in 
the L vs T diagram of Blain et al. (2004). We find that SSG 1 
is warmer than an equivalent luminosity SMG from the Chap- 
man et al. (2003; 2005) high redshift sample, and cooler 
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than the IRAS moderate redshift counterparts of Stanford et 
al. (2000). Instead, it is similar to the intermediate redshift 
ISO sample of GaiTett (2002). 

The notion of SSG 1 being an ISO ULIRG analog is further 
reinforced by the best-fitting 60-1 00 /im colors from STAR- 
DUST2, of which the closest template is FN 1-40, an ISO 
cold (Tdust=25.7K) ULIRG at z=0.45 (Chapman et al. 2002). 
The 850 /im flux (^3 mJy) predicts SSG 1 to be fainter than 
the majority of SCUBA-selected SMGs (Sgso > 5mJy, due to 
survey and instrument limitations). We conclude that SSG 1 
is more like /50-selected ULIRGs rather than the SCUBA- 
selected SMGs. 

6. CONCLUSION 

We report the detection of the first short-wavelength 
submillimeter-selected galaxy (SSG 1) and present optical, 
NIR and IR photometry which we use to constrain the red- 
shift, luminosity and dust temperature of this source. The 
photometric redshift estimators are all in agreement, produc- 
ing a best-fitting redshift of '-^l. The dust temperature of 
30.3±4.5K and luminosity of 1.0 x IO'^Lq make SSG 1 an 
analog of intermediate redshift /5(9-selected ULIRGs rather 
than high redshift SCUBA-selected SMGs. 

If there exists a population of galaxies with properties simi- 
lar to SSG 1, in a redshift space between ISO and SCUBA, 
it may be argued that observations at 350 /im, in conjunc- 
tion with the current ISO, Spitzer and SCUBA samples, could 
describe the IR-luminous population from 0.3 < z <4.0 (al- 
though the SCUBA-selected galaxies will be more luminous). 
Not only would this population bridge the gap in redshift 
space, but also the gap in "dust temperature space" (and con- 
sequently "dust mass space"). SHARC II thus complements 
SCUBA in revealing warmer SMGs that might be missed oth- 
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erwise. Furthermore, finding the lower luminosity SMGs is 
necessary for an accurate understanding of the submillimeter 
luminosity function. 

From our understanding of the first 350 /^m-selected object, 
SSG 1, we infer that observations at 350 /im, far from being 
superfluous to SCUBA, will be of great use in determining the 
nature and evolution of the luminous infrared galaxy popula- 
tion. 
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